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METHODS AND APPARATUS FOR TRUNCATION 
COMPENSATION 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims the benefit of U.S. provisional 
application number 60/416,072 filed October 4, 2002. 

BACKGROUND OF THE INVENTION 

[0002] This invention relates generally to methods and apparatus for 
computed tomographic (CT) image reconstruction, and more particularly to methods 
and apparatus for a truncation compensation scheme. 

[0003] Under some scanning conditions, portions of a patient may 
extend beyond the region measured by a detector, which may lead to image artifacts 
and an incomplete representation of the imaged object. Some known methods have 
been published that address the artifact reduction but not the imaging of the portion of 
the patient that is outside the field of view (FOV). However, it is desirable to image 
the portion of the patient that extends beyond the FOV. 

BRIEF DESCRIPTION OF THE INVENTION 

[0004] In one aspect, a method is provided. The method includes 
augmenting partially sampled field of view data using fully sampled field of view 
data. 

[0005] In another aspect, an imaging apparatus is provided. The 
imaging apparatus includes a radiation source, a detector responsive to radiation 
positioned to receive radiation emitted from the source, and a computer operationally 
coupled to the radiation source and the detector. The computer is configured to 
receive data from a Computed Tomography (CT) scan of an object, where the data 
includes fully sampled field of view data and partially sampled field of view data. 
The computer is further configured to augment the received partially sampled field of 
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view data using the fully sampled field of view data, and reconstruct an image of the 
object using the fully sampled field of view data and the augmented partially sampled 
field of view data. 

[0006] In a further aspect, a computer readable medium is provided 
which is encoded with a program configured to instruct a computer to augment 
partially sampled field of view data, and reconstruct an image using the fully sampled 
field of view data and the augmented partially sampled field of view data. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0007] Figure 1 is a pictorial view of a CT imaging system 

embodiment. 

[0008] Figure 2 is a block schematic diagram of the system 
illustrated in Figure 1 . 

[0009] Figure 3 illustrates truncation artifacts. 

[0010] Figure 4 is a graph showing a total attenuation integrated over 
all channels as a function of projection angle for a chest phantom. 

[001 1] Figure 5 is an illustration of truncation in a clinical setting. 

[0012] Figure 6 is a graph illustrating the impact of truncated 
projection on a total attenuation. 

[0013] Figure 7 is an illustration of slopes and boundaries estimation. 

[0014] Figure 8 is an illustration of fitted water cylinder for truncated 

projection. 

[0015] Figure 9 is an illustration of a projection extension scaled by 
the expected total attenuation. 

[001 6] Figure 10 illustrates a plurality of images. 
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DETAILED DESCRIPTION OF THE INVENTION 

[0017] There is herein provided truncation compensation methods 
and apparatus for Extended Field-of-View in Computed Tomography Systems. As 
explained in greater detail below, in one aspect, a method is at least partially based on 
a property that for parallel sampling geometry, the total amount of attenuation 
integrated over all channels is independent of the projection angle. The apparatus and 
methods are illustrated with reference to the figures wherein similar numbers indicate 
the same elements in all figures. Such figures are intended to be illustrative rather 
than limiting and are included herewith to facilitate explanation of an exemplary 
embodiment of the apparatus and methods of the invention. 

[0018] In some known CT imaging system configurations, a radiation 
source projects a fan-shaped beam which is collimated to lie within an X-Y plane of a 
Cartesian coordinate system and generally referred to as an "imaging plane". The 
radiation beam passes through an object being imaged, such as a patient. The beam, 
after being attenuated by the object, impinges upon an array of radiation detectors. 
The intensity of the attenuated radiation beam received at the detector array is 
dependent upon the attenuation of a radiation beam by the object. Each detector 
element of the array produces a separate electrical signal that is a measurement of the 
beam attenuation at the detector location. The attenuation measurements from all the 
detectors are acquired separately to produce a transmission profile. 

[0019] In third generation CT systems, the radiation source and the 
detector array are rotated with a gantry within the imaging plane and around the object 
to be imaged such that an angle at which the radiation beam intersects the object 
constantly changes. A group of radiation attenuation measurements, i.e., projection 
data, from the detector array at one gantry angle is referred to as a "view". A "scan" of 
the object includes a set of views made at different gantry angles, or view angles, 
during one revolution of the radiation source and detector. 

[0020] In an axial scan, the projection data is processed to 
reconstruct an image that corresponds to a two dimensional slice taken through the 
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object. One method for reconstructing an image from a set of projection data is 
referred to in the art as the filtered back projection technique. This process converts 
the attenuation measurements from a scan into integers called M CT numbers" or 
"Hounsfield units", which are used to control the brightness of a corresponding pixel 
on a display device. 

[0021] To reduce the total scan time, a "helical" scan may be 
performed. To perform a "helical" scan, the patient is moved while the data for the 
prescribed number of slices is acquired. Such a system generates a single helix from a 
fan beam helical scan. The helix mapped out by the fan beam yields projection data 
from which images in each prescribed slice may be reconstructed. 

[0022] As used herein, an element or step recited in the singular and 
preceded with the word "a" or "an" should be understood as not excluding plural said 
elements or steps, unless such exclusion is explicitly recited. Furthermore, references 
to "one embodiment" of the present invention are not intended to be interpreted as 
excluding the existence of additional embodiments that also incorporate the recited 
features. 

[0023] Also as used herein, the phrase "reconstructing an image" is 
not intended to exclude embodiments of the present invention in which data 
representing an image is generated but a viewable image is not. Therefore, as used 
herein the term, "image," broadly refers to both viewable images and data 
representing a viewable image. However, many embodiments generate (or are 
configured to generate) at least one viewable image. 

[0024] Figure 1 is a pictorial view of a CT imaging system 10. 
Figure 2 is a block schematic diagram of system 10 illustrated in Figure 1. In the 
exemplary embodiment, a computed tomography (CT) imaging system 10, is shown 
as including a gantry 12 representative of a "third generation" CT imaging system. 
Gantry 12 has a radiation source 14 that projects a cone beam 16 of X-rays toward a 
detector array 18 on the opposite side of gantry 12. 
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[0025] Detector array 18 is formed by a plurality of detector rows 
(not shown) including a plurality of detector elements 20 which together sense the 
projected X-ray beams that pass through an object, such as a medical patient 22. Each 
detector element 20 produces an electrical signal that represents the intensity of an 
impinging radiation beam and hence the attenuation of the beam as it passes through 
object or patient 22. An imaging system 10 having a multislice detector 18 is capable 
of providing a plurality of images representative of a volume of object 22. Each 
image of the plurality of images corresponds to a separate "slice" of the volume. The 
"thickness" or aperture of the slice is dependent upon the thickness of the detector 
rows. 

[0026] During a scan to acquire radiation projection data, gantry 12 
and the components mounted thereon rotate about a center of rotation 24. Figure 2 
shows only a single row of detector elements 20 (i.e., a detector row). However, 
multislice detector array 18 includes a plurality of parallel detector rows of detector 
elements 20 such that projection data corresponding to a plurality of quasi-parallel or 
parallel slices can be acquired simultaneously during a scan. 

[0027] Rotation of gantry 12 and the operation of radiation source 14 
are governed by a control mechanism 26 of CT system 10. Control mechanism 26 
includes a radiation controller 28 that provides power and timing signals to radiation 
source 14 and a gantry motor controller 30 that controls the rotational speed and 
position of gantry 12. A data acquisition system (DAS) 32 in control mechanism 26 
samples analog data from detector elements 20 and converts the data to digital signals 
for subsequent processing. An image reconstructor 34 receives sampled and digitized 
radiation data from DAS 32 and performs high-speed image reconstruction. The 
reconstructed image is applied as an input to a computer 36 which stores the image in 
a mass storage device 38. 

[0028] Computer 36 also receives commands and scanning 
parameters from an operator via console 40 that has a keyboard. An associated 
cathode ray tube display 42 allows the operator to observe the reconstructed image 
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and other data from computer 36. The operator supplied commands and parameters 
are used by computer 36 to provide control signals and information to DAS 32, 
radiation controller 28 and gantry motor controller 30. In addition, computer 36 
operates a table motor controller 44 which controls a motorized table 46 to position 
patient 22 in gantry 12. Particularly, table 46 moves portions of patient 22 through 
gantry opening 48. 

[0029] In one embodiment, computer 36 includes a device 50, for 
example, a floppy disk drive or CD-ROM drive, for reading instructions and/or data 
from a computer-readable medium 52, such as a floppy disk or CD-ROM. In another 
embodiment, computer 36 executes instructions stored in firmware (not shown). 
Generally, a processor in at least one of DAS 32, reconstructor 34, and computer 36 
shown in Figure 2 is programmed to execute the processes described below. Of 
course, the method is not limited to practice in CT system 10 and can be utilized in 
connection with many other types and variations of imaging systems. In one 
embodiment, Computer 36 is programmed to perform functions described herein, 
accordingly, as used herein, the term computer is not limited to just those integrated 
circuits referred to in the art as computers, but broadly refers to computers, processors, 
microcontrollers, microcomputers, programmable logic controllers, application 
specific integrated circuits, and other programmable circuits. Although the herein 
described methods are described in a medical setting, it is contemplated that the 
benefits of the invention accrue to non-medical imaging systems such as those 
systems typically employed in an industrial setting or a transportation setting, such as, 
for example, but not limited to, a baggage scanning CT system for an airport or other 
transportation center. 

[0030] Under some scanning conditions, portions of patient 22 
extend beyond the region measured by detector 18, which may lead to image artifacts 
and an incomplete representation of the imaged object. The X-ray tube and detector 
18 are rigidly held together on a frame which rotates around patient aperture 72. 
During the course of rotation, measurements are continuously made within a "fully 
sampled field of view" 60. The attenuation of x-rays traversing any areas of object 22 
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positioned outside fully sampled field of view 76 is measured at a limited range of 
rotation angles, and this region is referred to as the "partially sampled field of view" 
region. In other words, portions that are within fully sampled field of view 76 are 
positioned within fan 16 such that measurements are obtainable at all gantry angles, 
and the data collected is defined as fully sampled field of view data. Some portions, 
however, are within fan 16 at some angles but are outside fan 16 at other angles, and 
the data collected regarding these portions is defined as partially sampled field of view 
data. 

[0031] Some known methods have been published that address the 
artifact reduction but not the imaging of the portion of the patient that is outside the 
field of view (FOV). However, it is desirable to image the portion of the patient that 
extends beyond the FOV. This is useful in many fields including Oncology, Spin 
Angiography, Fused imaging systems, and In Economy CT Scanners. The current 
hardware of known multi-slice CT scanner limits the reconstruction field of view 
(FOV) to about 50 centimeters (cm.). Although this is sufficient for most clinical 
applications, it is desirable to expand the FOV to image objects outside this FOV. 
This may have particular advantages for applications such as oncology or CT/PET. 
For oncology applications, a larger FOV is desired. This is mainly due to the fact that 
for radiation treatment planning, the limbs of the patient are often positioned outside 
the scan FOV for better tumor positioning. The known CT reconstruction algorithms 
ignore the truncated projections and produce images with severe artifacts. These 
artifacts may affect an accurate estimation of the attenuation path for treatment 
planning. One phantom example is shown in Figure 3. For Fused Imaging systems 
such as CT/PET (Computed Tomography / Positron Emission Tomography) the FOV 
of the PET system may not match the existing CT design. It is desirable to have a 
consistent FOV between the CT and the other imaging system, CT/PET, CT/NUC 
(CT/Nuclear), or CT/MR (CT/Magnetic Resonence). This correction can be used to 
adjust the FOV to match. For PET this enables better attenuation correction. Herein 
described is an algorithmic approach to increase the reconstruction FOV beyond the 
FOV limited by the detector hardware. This correction algorithm can be applied to 
various reconstruction algorithms including, but not limited to, Full scan, 
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Halfscan/segment, helical scan, and the cardiac sector based algorithms. Additionally, 
system 10 is configured to employ the herein described algorithms. 

[0032] Figure 4 shows the total amount of attenuation plotted as a 
function of projection angle with a parallel beam sampling geometry for a chest 
phantom scan. Note that the curve is nearly a horizontal line. This property does not 
exist, however, for the fan beam sampling geometry. When the scanned object is 
outside the scan field-of-view (FOV), this property is no longer valid. The amount of 
deficit is equal to the portion of the object that is outside the projection FOV. In 
nearly all clinical cases, the projection truncation occurs only in a portion of the 
projection angles, as illustrated in Figure 5. In this example, the projection taken at a 
3 o'clock position is free of truncation and the projection taken at a 12 o'clock is 
severely truncated. Therefore, one can rely on the un-truncated projections (i.e., 
positions around 3 o'clock in Figure 5, for example) to estimate the amount of 
truncation for the truncated views (e.g., positions around 12 o'clock for the example in 
Figure 5). One early step in the correction process is to perform software fan-beam to 
parallel-beam rebinning on the pre-processed projections. In one embodiment, that 
early step is the first step. This process is well known in the art and does not require 
special data collection. Once the rebinning is completed, the projections are 
integrated over all detector channels to obtain the total attenuation curve, as shown in 
Figure 6. Note that the dips in the total attenuation curve correspond to views with 
truncation. The flat portion of the curve corresponds to the views in which no object 
truncation occurs. Once the total amount of the object outside the FOV is estimated, 
the next step is to estimate the distribution of the missing projection. To achieve this 
objective, in one embodiment, one first calculates the boundary reading, pi and p r as 
shown below in equation 1, in the truncated projection, as shown in Figure 7. To 
reduce noise, in one embodiment, the average of m number of samples is used. m=3 
has empirically been shown to be useful in reducing noise. In other embodiments, m 
is greater than 1 and less than 5. 
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[0033] Equation 1: pi= — £ and J* = ~~X p(N-i y k) 

m i=\ *n ;=1 

where N is the number of detector channels, and k is the projection view 

index. In addition, in one embodiment, the slopes, Si and Sr near both ends are also 

estimated. The slope estimation is performed by fitting n samples near the ends with a 

first order polynomial. n=5 has empirically shown to be useful. In one embodiment, 

n is greater than 2 and less than 8. In another embodiment, n is greater than 3 and less 

than 7. 

[0034] To further improve the reliability of the estimation, 
projections are used which are acquired from neighboring detector rows. Since 
human anatomy typically does not change quickly over a small distance (a few 
millimeters), the boundary samples and the slopes estimated from the neighboring 
rows do not typically vary significantly. Therefore, the estimated parameters (p ]9 p r , si, 
and Sr) can be the weighted average of the values calculated from several detector 
rows. Based on the boundary and slope information, one estimates a location and the 
size of a cylindrical water object that can be best fitted to the truncated projection. If 
we denote the attenuation coefficient of water /-iw, the radius of the cylinder R, and the 
distance from the cylinder center X, the projection value, p(x), and slope, p'(x), can be 
described by the following equation. 

[0035] Equation 2: p(x) = 2/j»Jr 2 - x 2 and p'(x) = " 2 ^ x 

VR 2 -x 2 

[0036] Since we calculate both p(x) and p'(x) at the truncated 
projection boundaries, the goal is to estimate R and x so that one obtains the size and 
location of the cylinder that is to be appended to the missing projection. The formula 
to estimate these parameters can be described by the following equations: 

[0037] Equation 3: xi = " ( f andRi = J^~t+ xi 2 
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[0038] Equation 4: Xr = A j and R, = J-^T + Xr ' 

[0039] The variables represent the estimated location and the size of 
the cylindrical objects that need to be extended from the truncated object. Once these 
parameters are determined, the expended projections can be calculated using equation 
(2). The process is illustrated in Figure 8. 

[0040] In this example, a cylindrical water phantom was used for 
simplicity. In reality, other object shapes, such as an elliptical cylinder, can also be 
used to increase the flexibility. If a priori information is available on the 
characteristics of the scanned object, the information can of course be used in the 
shape and material selection of the appending object. Iterative methods can be used to 
estimate the missing projection data. 

[0041] The estimated cylinders on both ends of the projection do not 
always recover the total amount of attenuation for the entire projection, since these 
objects are determined solely from the slope and boundary samples. None of the 
information derived from the total attenuation curve (Fig. 6) is used. To ensure proper 
compensation for the total attenuation loss, the attenuation distribution of the left side, 
Ti, verses the right side, T r , is determined based on the magnitude of pi and p r . 

[0042] Ti = PlT and 7V= — 

Pl+ pr pi\ pr 

[0043] where T is the total amount of loss in attenuation determined 
from Fig. 6. Additionally, if the amount of attenuation under the extended curve is 
insufficient to make up for the attenuation loss, the estimated projection is stretched to 
fulfill the attenuation deficits, as illustrated in Figure 9. On the other hand, if the 
amount of attenuation under the extended curve exceeds the attenuation loss, the 
estimated projection is compressed in a similar fashion. In one embodiment, the 
calculation process is as follows. One first calculates the ratio of the expected total 
attenuation (shown in Eq. (5)) over the area under the extended projection curve 
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(shown by the shaded region in Fig. 9). If the ratio is larger than unity, the x-axis is 
scaled by the ratio so that the initial estimated projection (shown by the dotted line in 
Fig. 9) is further expanded (shown by the solid thick line in Fig. 9). Similarly, if the 
ratio is significantly smaller than unity, the expanded projection can be compressed in 
x. 

[0044] Figure 10 shows an example of the reconstructed phantom 
images without and with correction. A shoulder phantom was scanned in an axial 
scan mode with a 4 x 1.25 mm detector configuration. A 15cm plastic phantom was 
attached to the shoulder phantom in such a way that the edge of the plastic phantom is 
near the boundary of the 65cm FOV. The truncated object is nearly completely 
recovered. Note that Fig. 10(a) was reconstructed with a 50cm FOV without 
truncation correction (current product limit) and Fig. 10(b) was reconstructed with 
65cm FOV with the herein described correction. For reference, the phantom that is 
partially truncated is in shown in Fig. 10(c). 

[0045] Although the above described system and methods uses only 
the conservation of total attenuation, the magnitude, and the slope of the boundary 
samples to estimate the missing projection distribution, additional information can 
also be used for the estimation. For example, one could use the Helgason-Ludwig 
Condition (HL condition) for tomography to further refine the above technique. 
Additionally, different thresholds can be placed to ensure that the algorithm functions 
properly under erroneous measurement conditions. For example, one could set the 
upper and lower limits on the stretch ratio described in Fig. 9 to prevent the condition 
of increased error due to unreliable measurement. In addition, the slope calculation of 
si and Sr can be set such that it falls within a reasonable range. If the characteristic of 
the material of the scanned object is known to be significantly different from water, 
one can also use the attenuation coefficients of the known material (instead of water) 
to perform size and location calculations shown in Eqs. (3) and (4). 

[0046] Because the Interpolated data does not have the same image 
quality as data within the fully sampled FOV it may be useful to tag the image where 
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the FOV becomes extrapolated. In one embodiment, a delineation is provided in a 
reconstructed image between areas representative of the fully sampled field of view 
data and the partially sampled field of view data. Figure 10(d) illustrates the boundary 
is tagged by a dotted line. This could also be done with a color code or shift in the CT 
Number. Because the tag may impact the ability to view the image data an easy way 
is provided to turn on and off the tagging. A user of system 10 is allowed to turn on 
or off the tagging. 

[0047] While the invention has been described in terms of various 
specific embodiments, those skilled in the art will recognize that the invention can be 
practiced with modification within the spirit and scope of the claims. 
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